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ABSTRACT 

Most of the mass density in the Universe — and in the halo of our own 
galaxy — exists in the form of dark matter. Overall, the contribution of lu- 
minous matter (in stars) to the mass density of the Universe is less than 
1%; primordial nucleosynthesis indicates that baryons contribute between 
1% and 10% of the critical density (O.Ol/i" 2 < Q B < 0.02/i~ 2 ; h = the Hub- 
ble constant in units of lOOkms -1 Mpc" 1 ); and other evidence indicates that 
the total mass density is at least 10% of critical density, and likely much 
greater. If the universal density is as low as 10% of the critical density there 
may be but one kind of dark matter. More likely, the universal density is 
greater than 10%, and there are two kinds of dark matter, and thus two dark 
matter problems: In what form does the baryonic dark matter exist? and 
In what form does the nonbaryonic dark matter exist? The MACHO and 
EROS collaborations have presented evidence for the microlensing of stars in 
the LMC by 10~ M dark objects in the halo of our own galaxy and may 
well have solved one of the dark matter puzzles by identifying the form of 
the baryonic dark matter. It is too early to make precise statements about 
the fraction of the mass density in the halo of our galaxy contributed by 
lensing objects (= f m ), though the EROS/MACHO data suggest that f m is 
probably 0.1 or larger. Taking our galaxy to be typical and taking account 
a fraction f m of the mass in the portion of the halo that contributes most 
significantly to microlensing (within 20kpc of the galactic center), I estimate 
that lensing objects contribute a fraction 0.008 f m /h of the critical density, 
and clearly cannot account for the bulk of the dark matter if Qq 3> 0.1 (even 
discarding the nucleosynthesis bound). (If the distribution of lensing objects 
extends throughout the galactic halo, out to about 100 kpc, their contribution 
is larger by a factor of five). If halo lensing objects contribute around 10% 



of the mass density in the halo, the EROS/MACHO results provide further 
evidence for a nonbaryonic component to the halo dark matter, consistent 
with the simplest prediction for the ratio of cold dark matter (e.g., axions 
or neutralinos) to baryons in the halo. In any case, the EROS/MACHO 
results in no way lessen the strong case that now exists for the mass density 
of the Universe being significantly larger than that which baryons can con- 
tribute; nor do they affect significantly the prospects for directly detecting 
particle dark matter in our vicinity. Discovering the nature of the baryonic 
dark matter should provide still further impetus for solving the more weighty 
dark matter problem, the nature of the nonbaryonic dark matter. 



1 Introduction 



1.1 The events 

Two collaborations searching for massive, dark objects in the halo of our 
galaxy through gravitational microlensing reported candidate events this 
week. The French EROS collaboration has monitored the brightnesses of 
3 million stars in the Large Magellanic Cloud (LMC) over a three-year pe- 
riod and reported two events: a star that brightened by a factor of about 2.5 
over a characteristic time interval of about 54 days and a star that brightened 
by a factor of about 3.3 over a time interval of about 60 days [[[]]. (EROS 
measure the brightnesses of stars on 5° x 5° Schmidt plates taken with two 
filters, bleu and rouge, with a machine known as "MAMA.") The American- 
Australian MACHO collaboration has monitored the brightnesses of about 
1.8 million stars in the LMC for one year and reported one event: a star that 
brightened by about a factor of about 6.8 over a characteristic time interval 
of about 34 days [|2| .[] (MACHO use two CCD cameras and a dichroic beam- 
splitter to measure the brightnesses of stars in two colors, red and blue, and 
observed each star about 250 times during the year period.) 

1.2 Microlensing 

The basic idea of detecting massive, nonluminous objects in the halo by mi- 
crolensing was suggested by Paczynski || ; the discussion that follows is based 
upon the very nice paper by Griest |3Q. The actual deflection of light from a 
star in the LMC (distance 50 kpc) by an object in our halo is far too small to 
measure, of order 5(f) ~ 2 x lO _4 (m/M ) 1//2 arcseconds, where m is the mass 
of the halo object and the impact parameter (distance between the lensing 
object and the line of sight to the star) is taken to be the Einstein radius 
(see below). However, due to gravitational focussing the total brightness of 
the two unresolved images is greater than that of the unlensed star, by a 
factor equal to 1.34 for an impact parameter equal to the Einstein radius 
(and larger for smaller impact parameters). Since the a priori brightness of 
a given star is not known, Paczynski's idea was to look for time variation of 
the brightness of the star due to the motion of the lensing object across the 
line of sight. 

1 Thc definitions of event duration used by MACHO and EROS differ by about a factor 
of two; I have used the MACHO definition for all three events. 



1 



The brightness amplification A depends upon the impact parameter I in 
units of the Einstein radius Re- 



A ^ + 2 

u {u2 + 4)1/2' 

where / = uRe and 
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Here x is the distance to the lensing object and L is the distance to the star. 
A plot of A[u) is shown in Fig. 1 (Fig. 2 of Ref. @]); large A corresponds to 
small u. 

Because the lensing object is moving (the rms halo velocity is around 
300 km s -1 ) the impact parameter changes with time, and hence the amplifi- 
cation does too (in principle, one also must take into account the velocity of 
the star and the solar system, though it is a small correction [|J]). The time 
profile of the amplification, or light curve A(t), depends upon u(t) alone. 

To be specific, the impact parameter u is given by 

u{t) 2 = < in + (v T t/R E ) 2 , (3) 

where w m i n is the minimum value of the impact parameter, vt is the velocity 
of the lensing transverse to the light of sight, and the epoch of maximum am- 
plification defines time zero. Note that the light curves are a two-parameter 
family of curves. By the measuring the light curve for a microlensing event 
one can infer u min [= u(A max )] and Re/vt- Recall, Re ~ GmL/c 2 , which 

implies that the event duration t ~ V GmL / cvt ~ 100 da-y/m/M (more 
later). 

EROS define the duration of the event to be the characteristic time 
^eros = Re/vt, which corresponds to impact parameter u(±t ERO s) = v/l + Wmin! 
MACHO define the duration of the event to be the total time the amplifica- 
tion is above threshold for detection of amplification, £macho = 2^Ju^ — v%~Re/vt- 
For maximum amplification A max 3> 1 and a threshold amplification At ~ 
1.34 (ut ~ 1), which seems to apply to the data of both collaborations, 

^MACHO — 2tEROS- 

Now the formulas relating event rate to the abundance of halo objects; 
to begin, assume that the lensing objects have mass m, are distributed like 



(1) 
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the halo material, and contribute a fraction f m of the halo density, 



Prn 



f, 




(4) 



Here r = distance from the center of galaxy, a = 5 ±5 kpc is the core radius of 
halo, r — 8.5 kpc is the distance of the Solar System from the galactic center, 
and piocai — 5 x 10 _25 gcm -3 is the local halo density. (The mass density in 
the halo is determined by measurements of the rotation curve of our galaxy 
and the amount and distribution of light in the galaxy; the full extent of 
the halo — and its total mass — are not known, though there is evidence that 
the halo extends at least as far as the LMC and probably out to 100 kpc 
JJ.) Further, let us assume that the threshold for detection is amplification 
At, corresponding to impact parameter ut, and that the efficiency for the 
detection of a microlensing event above the amplification threshold is e. (The 
latter is clearly an oversimplification as e will be a function of A, the duration 
of the event, and perhaps other things.) 

The probability that a given star in the LMC is being microlensed with 
amplification greater than At by a halo object is referred to as the optical 
depth for microlensing, and was calculated by Griest [|J, cf. Eq. (4), 



where the halo is assumed to extend out at least as far as the LMC, the core 
radius is taken to be 5 kpc, ut is the impact parameter corresponding to At 
(ut = 1 for At = 1.34). The value of r is not very sensitive to the first two 
assumptions. 

Another interesting quantity is r(y), the probability that a star in the 
LMC is microlensed by a halo object between the galactic center and dis- 
tance y from it; r(y) is shown in Fig. 2. (Said another way, r(y) is the 
probability for microlensing under the assumption that the halo only ex- 
tends out to radius y.) Figure 2 illustrates that most of the optical depth for 
microlensing is due to halo objects within 10—20 kpc or so of the center of the 
galaxy.Q That is, microlensing only probes the inner portion of the halo, and 
further, not much of the lensing probability is due to objects near the solar 

2 The reason is simple: the Einstein radius, which the sets the radius of the tube within 
which halo objects lead to amplification above threshold, is proportional to the product 
of the distance to the lens and the distance from the lens to the LMC, and further, 
dr/dy oc p m /m, which falls off as r~ 2 . 



T ~ 4 X 10 7 f m U T ; 



(5) 
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neighborhood. This is an important fact to keep in mind when interpreting 
the EROS/MACHO events. 

What makes microlensing observable is the time variation of the ampli- 
fication. So even more important than the lensing probability is the rate of 
microlensing events (roughly r divided by the typical event duration); inte- 
grating over the distribution of velocities and positions of halo lensing objects 
Griest obtains the microlensing rate (per star observed) with amplification 
greater than At, cf. Eqs. (7-14), 



Note, T, unlike r, depends upon the mass of the halo objects: for fixed f m , 
the smaller the mass, the higher the event rate. 

The duration and amplification of a microlensing event depends upon 
the impact parameter and the velocity, mass, and distance of the halo mi- 
crolens. However, the velocity and distance of the halo lensing object can 
only be specified statistically, and thus the mass of the lensing objects can 
only be inferred from a measured light curve in a statistical way. Griest has 
constructed the likelihood function for the lensing mass, given the event du- 
ration and threshold for detection. The distribution is very broad, about a 
order of magnitude in mass at full- width, half maximum (Fig. 9 in Ref. H). 
The peak of that likelihood function occurs at a mass 



Note that this mass depends upon the threshold of the experiment; for EROS 
and MACHO u T ~ 1. 

The distribution of maximum amplification of microlensing events is basi- 
cally geometric: The fraction of events with maximum amplification greater 
than A corresponds to events with minimum impact parameter less than 
Wmin = u(A), which is simply equal to u(A)/ut- This fraction e(A) can also 
be expressed in terms of A and At, 



and is show in Fig. 3; cf. Fig. 10 of Ref. Q. 

To conclude this section; microlensing has many clear signatures: it is 
achromatic; the light curve is defined uniquely in terms of two parameters, 




(6) 




(7) 




(8) 
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maximum amplification and duration; the distribution of event amplifications 
is predicted; microlensing shows no preference for the type of lensed star; and 
so on. The most significant background for microlensing searches are variable 
stars; already, MACHO has compiled the largest catalogue of variable stars in 
the LMC. Because of its many signatures, microlensing should be a relatively 
easy hypothesis to test; indeed, EROS and MACHO have already mastered 
the difficult tasks of monitoring the brightnesses of millions of stars and 
discriminating against known types of variable stars. In the next year or 
so both collaborations should increase their data sets by a factor of ten or 
so and thus should clarify the few questions that have been raised about 
their events (large x 2 f° r two of the three events and large amplification for 
the MACHO event despite their low threshold, cf. Fig. 3). The case for 
(or possibly against) the existence of halo lensing objects should be firmly 
established soon. 

2 Overinterpreting the Events 

What can one learn from the very preliminary MACHO and EROS results? 
Definitely not as much as one would like! But let's risk going out on a limb 
with some shaky interpretations. 

First, the number of events seen. Based upon their stated exposures the 
number of expected events is 



where of course the detection efficiencies and thresholds for the two exper- 
iments are not necessarily the same. In order to infer the fraction of the 
mass density of the halo in lensing objects, f m , one must know not only the 
efficiencies, but also something about the mass of the lensing objects (see 



MACHO has made no statement about f m . From the duration of their 
event they infer a most likely mass of about O.1M ; from their paper and 

3 What one wants to do, and what both collaborations are certainly doing, is construct- 
ing the likelihood function for their events as a function of all the relevant parameters 
(lens mass or distribution of lens masses, spatial density of lenses, and so on) and taking 
account of acceptance, etc. From this, the most likely values for m and f m and estimates 
for their variances can be inferred. 



-^MACHO 

-Weros 



9.4/ m e MT /(m/O.lM ) 1 / 2 
47/ m ^ T /(m/O.lM ) 1 / 2 



(9) 
(10) 



below) .0 
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conversations with collaboration members I conclude that both ut and e are 
of order unity. Climbing out to the very end of the limb, this suggests that 
f m is about 0.1 or larger (with great uncertainty). 

EROS state in their paper that "the number of events is consistent with 
the hypothesis that the halo of our galaxy is made essentially of dark objects 
in the mass range of a few 10~ 2 to a few M ." Taking m = 0.1M Q , this 
implies an efficiency cut ~ 0.04. My interpretation of a talk given by EROS 
member Luciano Moscoso is that, based upon Monte Carlo simulations, the 
detection efficiency for events with w min < 1 is about 80% and for events 
with u m in > 1 it is much, much smaller. Their exposure seems to be much 
lower than 3 million stars times 3 years: only about half the stars were used 
and the time coverage of a typical star was about 600 days. Putting this 
together, I guess that ut ~ 1 and e ~ 0.1, suggesting that f m ~ 0.2 or so. 

Next, event duration; as mentioned earlier, the likelihood function com- 
puted by Griest peaks at a mass of 1.5M (t/lOO da) 2 /i^. For MACHO, this 
implies a lensing mass of about 0.1M o , which is the central value stated in 
their paper. For EROS, the event durations are about twice as long; this 
implies a most likely mass of about O.4M . While the estimates for the 
lensing mass seem to be different, recall that the likelihood function is an 
order of magnitude in mass at FWHM and that the detection efficiencies will 
certainly depend upon the event duration and can influence the likelihood 
function. 

As noted above, both the event rate and event duration are needed to infer 
the fraction of critical density contributed by halo lensing objects. Moreover, 
some a priori information may help in this regard. For example, stars made of 
hydrogen and helium more massive than about O.1M are nuclear burning; 
previous unsuccessful searches for low-mass stars (M dwarfs) place severe 
limits to f m for m > O.1M , perhaps even ruling out this possibility || [/]]. 
Dark halo objects heavier than around a solar mass would have to be neutron 
stars, black holes, or white dwarfs, otherwise they would be easily visible; 
however, this possibility too is severely constrained by their contribution to 
the heavy-element abundance in our galaxy 0. For these reasons it has 
been argued very convincingly that the mass of baryonic halo objects must 
be less than about O.1M ; note, this makes my previous shaky estimates for 
f m shaky upper limits since Y oc m -1 / 2 and f m oc m 1 / 2 . 

Let's go on and estimate the fraction of critical density contributed by 
halo lensing objects, fl m . The strategy is as follows: (1) construct the ratio 
of mass in halo lenses to light for the Milky Way, M m f L; (2) assume that 
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this ratio is universal; (3) Divide this ratio by the critical mass to light ratio 
to find Q m .f\ 

Much is known about the halo mass interior to about 20 kpc ||, and recall, 
most of the microlensing is due to objects within 20 kpc of the galactic center. 
The halo mass interior to 20 kpc is about 2.3 x 1O 11 M . The luminosity of 
the Milky Way galaxy is about 2.3 x 1O 1O L B0 . Assigning a fraction f m of the 
halo mass within 20 kpc of the galactic center to lensing objects, I find 



r Wmy^- (11) 

Q m ~ 8 x 10- 3 f m /h. (12) 

This is not much mass, though, to be sure, at least one of our assumptions 
was very conservative. Although microlensing searches are most sensitive to 
halo lensing objects within 20 kpc of the galactic center, it would be surprising 
if f m suddenly dropped to zero at this point. If we assume that halo lensing 
objects populate the halo at the same fractional abundance out to 100 kpc, 
the known extent of the halo, the previous estimate increases by a factor of 
five: 

Q m ~ OMfJh, (13) 

though still far from closure density. 

To conclude this discussion, the MACHO and EROS collaborations may 
well have solved one of the dark matter riddles, namely, what form the dark 
baryons take. If so, their searches were spectacularly successful. However, 
the weightier dark matter problem is the nature of the nonbaryonic dark 
matter, as evidence indicates that Qq exceeds 0.1, the maximum contribution 
permitted for baryons. (Nucleosynthesis bound aside, the contribution of 
halo lensing objects to the universal density cannot account for the total 
mass density if Qq is significantly greater than 0.1.) 



3 Implications for Particle Dark Matter 

There is mounting evidence that Qq is significantly greater than 0.1, the 
maximum contribution of baryons, and perhaps close to unity. The evidence 

4 Astronomers measure the mass density of the Universe by multiplying the luminosity 
density, C ~ 2 Ah x 10 8 £b© Mpc -3 , times a characteristic mass to light ratio; the mass to 
light ratio (in the Bt system) that corresponds to critical density is 12OO/iM /Lb . 
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that f2o is close to unity includes several studies relating the peculiar ve- 
locities of galaxies within a few 100 Mpc of the Milky Way (and the Milky 
Way itself) to the distribution of galaxies determined from red-shift sur- 
veys based upon the IRAS Point Source Catalogue^ and the fact that the 
only models for the formation of structure in the Universe that are consistent 
with measurements of Cosmic Background Radiation (CBR) anisotropy (e.g., 
COBE and experiments on smaller angular scales) and large-scale structure 
require nonbaryonic dark matter. A number of measurements of Qq, e.g., 
cluster virial masses and the infall of nearby galaxies into the Virgo Cluster, 
strongly suggest that Qo is greater than 0.1, though perhaps not as large as 
unity. And finally, there are theoretical reasons in favor of Qq = 1 and non- 
baryonic dark matter: a flat Universe is an unambiguous prediction of the 
inflationary paradigm, and there are several particles whose postulated exis- 
tence is motivated by compelling particle-physics considerations and whose 
relic abundance is close to the critical density (e.g., neutralino and axion) 

In a flat, critical Universe particle dark matter is obligatory. For the 
discussion that follows I assume that Qq = 1 and that Qb ~ 0.05 — 0.10. The 
universal ratio of exotic dark matter to baryonic dark matter is then 

^ = i^UlO-20. (14) 

Of greater interest for the interpretation of the MACHO/EROS events is the 
local ratio of particle dark matter to baryonic dark matter. In order to fully 
answer this question we must know how galaxies and other structures formed. 
In the simplest circumstance, where only gravity is important and the initial 
velocities of baryons and particle dark matter are similar (or negligible), the 
answer is dictated by the equivalence principle: Baryons and particle dark 
matter must follow the same trajectories so that the ratio remains constant 
and equal to 10 — 20. 

When nongravitational forces become important then the two forms of 
dark matter will certainly become differentiated as particle dark matter does 
not interact electromagnetically, the most important nongravitational force 
that operates in the Universe. In general, the interactions of baryonic matter 



5 The peculiar velocities of galaxies arise due to gravitational forces resulting from the 
inhomogeneous distribution of matter and depend upon the distribution of galaxies and 
the mean mass density; by determining the distribution of galaxies one can infer the mean 
mass density; see e.g., Refs. [0. 
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with itself (and the CBR) allow baryons to "cool" (i.e., dissipate their energy) 
and become more condensed, thereby decreasing the ratio of particle dark 
matter to baryonic dark matter. 

For hot dark matter, relic particles that move fast (e.g., 20 eV neutri- 
nos), the velocity distribution of the particle dark matter is very important: 
Neutrinos move too fast to become trapped in galaxies (at least until very 
recently) and so one expects the local ratio of hot dark matter to baryonic 
matter to be much, much lower than the universal ratio of 10 — 20. 

For cold dark matter (e.g., axions or neutralinos) , relic particles that move 
very slowly, the velocity distribution is unimportant and one would expect 
that in objects whose formation only involves gravity the local density would 
reflect the universal value of 10 — 20. There is no doubt that the formation 
of the disk of our galaxy involved dissipation (the disk-like structure traces 
to the dissipation of essentially all the energy not associated with angular 
momentum). While the halos of galaxies are not fully understood, they show 
no obvious signs of dissipation having been involved in their formation. Thus, 
the simplest assumption is that the ratio of cold dark matter to baryonic dark 
matter in the halo should reflect the universal value of 10 — 20.0 

The fraction of the halo mass density in lensing objects is crucial to un- 
derstanding the implications of MACHO and EROS for particle dark matter. 
Let us consider two possibilities consistent with the data at hand: f m ~ 0.1 
and f m ~ 0.5. 

Suppose f m turns out to be of order 0.1, as is suggested by a naive over- 
interpretation of the data. In this case the MACHO/EROS results are con- 
sistent with the simplest expectation for cold dark matter and inconsistent 
with hot dark matter. Moreover, the fact that 90% of the mass density in the 
halo remains unexplained would provide even further impetus for searching 
for cold particle dark matter in our own halo. 

On the other hand, suppose that two or three years from now EROS and 
MACHO convincingly establish that f m is 0.5. This is consistent with hot 
dark matter and inconsistent with the simplest expectation for cold dark 
matter. With regard to the latter I mention: 

1. It could be that the inner part of the halo (which is what is probed by 
MACHO and EROS) has a smaller ratio of cold dark matter to baryons 

6 The ratio of cold dark matter to baryonic dark matter expected for the mixed dark 
matter model (fts ~ 0.05, f2 C oid ~ 0.65, f2hot ~ 0.3) is only slightly smaller, around 7— 13. 
This is one of the models that best accounts for the observed large-scale structure. 
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because the baryons underwent some dissipation (after all they cannot 
have formed into astrophysical objects without dissipation). At present 
there are no compelling scenarios for the formation of these objects, so 
it is certainly possible that the baryons in the inner halo have undergone 
some dissipation. 

2. It could be that the universal ratio of cold dark matter to baryons is 
only a few. One of the models that best accounts for the observed large- 
scale structure is cold dark matter + cosmological constant (f^A = 0.8; 
^coid ~ 0.15; and Qb ~ 0.05). While it is certainly not the best 
motivated model, if the Hubble constant is close to 80kms _1 Mpc -1 , 
as many now believe ||, a cosmological constant will be mandatory to 
solve the cosmological age crisis. 

3. A halo ratio of cold dark matter to baryonic dark matter of order unity, 
rather than 10 — 20, has only a minor impact on experimental searches 
for the cold dark matter in our vicinity, as these experiments already 
take into account a factor of two uncertainty in the local density of halo 
material .[] 

4 Concluding Remarks 

Luminous matter (in stars) contributes less than 1% of the critical density, 
and primordial nucleosynthesis makes a strong case that ordinary matter 
contributes between 1% and 10% of the critical density, establishing the 
existence of the first dark matter puzzle, the form of the baryonic matter. 
There is compelling evidence that the mass density of the Universe is much 
greater than 10% of the critical density and perhaps close to the critical 
density, establishing the existence of a second dark matter puzzle, the form 
of the nonbaryonic dark matter. 

7 Given the uncertainties about the halo mass, even within 20kpc of the galactic center, 
and the insensitivity of microlensing to the density of lensing objects nearby, it seems 
unlikely that microlensing experiments could ever establish with confidence that lensing 
objects contribute more than about 50% or so of the local dark matter density. Moreover, 
because there is nothing to prevent cold dark matter from falling into the halo of our 
galaxy, even if baryonic matter has been concentrated by dissipation, it is difficult to 
imagine that the local density of cold dark matter could be significantly less than the local 
halo density. 
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The EROS/MACHO collaborations may well have solved the first dark 
matter problem by providing evidence for dark baryons in the form of as- 
trophysical mass objects in the halo of our own galaxy.^ However, their 
discovery, as important as it is, has little to say about the second dark mat- 
ter problem: The mass density contributed by halo lensing objects cannot 
significantly exceed 10% of the critical density (even discarding the constraint 
to Qb based upon primordial nucleosynthesis). 

A critical Universe is well motivated, supported by mounting evidence, 
and must involve particle dark matter. The detection of dark stars in the halo 
of our own galaxy changes none of this; nor does it lessen significantly the 
prospects for the direct detection of the cold dark matter in our vicinity (e.g., 
in laboratory experiments designed to detect halo axions or neutralinos) . The 
scientific stakes in testing the particle dark matter hypothesis are extremely 
high: identifying the primary constituent of the Universe, discovering a new 
particle of nature, and probing the earliest moments of the Universe. The 
pioneering searches for particle dark matter |TT[ must continue at full speed; 
if anything, solving part of the dark matter puzzle should provide still further 
impetus for solving the rest of the puzzle! 
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ments. This work was supported in part by the Department of Energy (at 
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5 Figure Captions 



Figure 1: Impact parameter u (in units of the Einstein radius, u = I/Re) 
function of the amplification. 

Figure 2: The optical depth for gravitational lensing due to halo objects 
within a distance y of the galactic center, r(y), as a function of y. Note that 
most of the optical depth for microlensing is due to objects within 10 — 20 kpc 
of the galactic center. Thus, microlensing experiments only probe this part 
of the halo and are insensitive to the extent of the halo as well as the local 
halo density. 

Figure 3: The fraction of microlensing events expected with amplification 
greater than A for amplification thresholds of A T = 1.1, 1.34, 1.5, 2, 3, and 5. 
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